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Abstract 

We study directed random graphs (random graphs whose edges are directed) as they 
evolve in discrete time by the addition of nodes and edges. For two distinct evolution 
strategies, one that forces the graph to a condition of near acyclicity at all times 
and another that allows the appearance of nontrivial directed cycles, we provide 
analytic and simulation results related to the distributions of degrees. Within the 
latter strategy, in particular, we investigate the appearance and behavior of the 
strong components that were our subject in the first part of this study. 
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1 Introduction 



In the first part of this study [1], we considered random digraphs on a fixed 
set of n nodes with Poisson-distributed in- and out-degrees. For a wide range 
of values for z (the mean in- or out-degree of a node), we investigated the 
behavior of the strong components of the digraph, which essentially are max- 
imal subgraphs whose nodes can all be reached from one another by following 
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the directions of the edges. We contributed analytic and simulation results to 
the special cases of cycle components and knots. More specifically, we demon- 
strated that cycle components are concentrated near z = 1 and necessarily 
small, encompassing a number of nodes proportional to In n, while knots tend 
to occur only sparsely as size-one components for very small values of z, then 
become practically absent as 2; is increased, then finally occur again as a sin- 
gle knot that encompasses nearly all n nodes for larger values of z. The latter 
transition is sharp and happens roughly at z = ln(2n). 

In this second part, our aim is to study the behavior of strong components 
when the digraph is no longer static, but rather evolves in time. This type of 
study was also pioneered by Erdos and Renyi [2], who considered undirected 
graphs on a fixed set of nodes that progressively become more and more in- 
terconnected by the random addition of edges between pairs of nodes, thus 
leading to a Poisson distribution of node degrees. Recently, though, in an at- 
tempt to model the networks that occur in many areas of interest, considerable 
effort has been directed towards studying evolution scenarios in which nodes 
and edges may enter and leave the network continually. The reader is referred 
to the surveys in [3] for details related to several areas. One of the surveys 
[4] highlights the mathematics of a variety of such models (e.g., [5]), many of 
which exhibit the power-law (as opposed to Poisson) distribution of degrees 
that is characteristic of networks that grow by some sort of preferential (as 
opposed to random) attachment of nodes to each other. 

Our focus is on the study of digraphs that evolve by the addition of new 
nodes, as well as new edges, in discrete time. Also, although we do touch the 
issue of preferential attachments briefly at one point, for the most part we 
follow [1] and concentrate on the case of random connections. Even though 
some of the models discussed in the literature (cf. [3]) can be claimed to 
include digraphs — at least inherently — we believe our emphasis on random 
connections in the evolution of digraphs covers new ground. 

The following notation is common to all the remaining sections. For t > an 
integer, we let Df denote the digraph at time t and Nt, with rit = \Nt\, its 
node set. Dq is assumed to have no nodes, so A^^q = and Uq = 0. Unlike the 
static case [1], it is now necessary to consider in- and out-degrees separately. 
For node i, we let d^{i) denote its in-degree at time t and df{i) its out-degree, 
with expected values Zf'{i) and z^{i), respectively. 

The paper is organized as follows. In Section 2 we consider evolution scenarios 
in which the deployment of edges disallows every nontrivial directed cycle. In 
Section 3 we modify the deployment rule in order to allow the appearance 
of arbitrary directed cycles, and with them nontrivial strong components. 
Conclusions are given in Section 4. 
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2 NeEirly acyclic evolution 



2.1 Analytic results 



In the first evolution scenario that we consider, at every time step t > a 
new node is added to the digraph, so we have = t. We assume that nodes 
are numbered consecutively from 1 as they enter the digraph, so node i is the 
node added at time t — i. When a node enters the digraph, a random number 
of edges is also added, all of them directed away from the newly added node. 
Except for the possibility of self-loops (an edge leading from a node to itself 
is called a self-loop), the digraph is at all times acyclic (i.e., has no directed 
cycles spanning more than one node) and all its strong components have one 
single node. 

For some fixed z, we start by considering the case in which an edge is deployed 
from node i to each of nodes 1, . . . ,i independently with probability z/i (pro- 
vided, of course, that z < i; for z > i, it must be assumed that i is connected 
to all of 1, . . . , i). In this case, for sufficiently large i and t > i, the out-degree 
df{i) is clearly Poisson-distributed with mean zf{i) — z. 

As for the in-degree df{i), its expected value zf{i) is, for sufficiently large i 
and t > i, given by 

zt{i) = j:"- = z{H,-H,.,\ (1) 

u=i 



where H^n = J2u=i "^^^ harmonic number, given for large m by 

Hjji = Inm -I- 7, with 7 denoting Euler's constant [6]. Thus, for i » 0, 

ztii) ^ In (ij. (2) 



In order to discover how df~{i) is distributed, let Pt{i. k) denote the probability 
that d^{i) = k for t > i and k > 0. li t — i + 1 < k, then clearly Pt{i, k) = 0. 
Otherwise, ioi k <^t — i and z ^ i,t, we have 



Pt{t,k) 



Pt{i,0) 
k\ 



(3) 



as we demonstrate in Appendix A. Pt{i, 0) can now be approximated from (3) 
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by summing Pt{i, k) over the appropriate range of k: 



t-i+l 1 . 

1 ^ PS, 0) E ^ k«l ^ Ptih 0)e^*"», 



k=0 



(4) 



hence Pt{i, 0) ~ e ^t^*^*^ and (3) can be re-written as 

m,k) 





k 







A;! 



(5) 



which is the Poisson distribution with mean zfii). 

Denoting the overall in-degree distribution by Pt{k), we see that it can be 
estimated by averaging (5) over all nodes, that is, 



^ i=l 



z\nt 



k\ 



4t / a;'=e-(^+V^)-rf^^ 
zk\ J 



(6) 



x=0 



where x = z\n(t/i) (hence di = —{t/ z)e~^^^dx). The integrand in (6) peaks at 
X = kz/{z + 1) with a finite value, so the integral can be extended to infinity, 
yielding 



CO 

Ptik) ~ ^ / x'^e-^^+^l'^^'dx 



(^ + 1) 



fe+1 ' 



(7) 



Let us now turn to the case in which node upon entering the digraph at time 
t — connects out to node j e {1, . . . , i} with probability proportional to how 
many incoming edges j already has (plus 1, to ensure nonzero probabilities to 
start with) while aiming at the same mean out-degree z. In this case, an edge 
is deployed from i to j with probability 



(8) 



with the understanding that df_i{i) — 0, provided z[l + df_i{j)] < t + 
J2u=i dtLiiu) (when this does not hold, then the deployment is assumed to 
take place with probability 1). While (8) trivially ensures that the expected 
value z^{i) of d'i~{i) is equal to z for all t that render the expression in (8) a 
legitimate probability, the underlying distribution of the out-degree df{i) is 
not so easily deduced, although by the results surveyed in [3] its average at 
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time t over all the t nodes already present in the digraph must be proportional 
to k~°' for some a > 0. 

The same holds for the average distribution Pt{k) of the in-degrees df{i). The 
expected value Zt~{i) of df(i), in turn, can be assessed as follows. By (8), the 
expected probability that node i connects out to node j at time t — i is 



Thus, for i > i, Zt'{i) obeys the recurrence 



z^{i) = zl_^{i) + z 



' l + ^tti(i) ' 



H-\\ 

z 



V ^(-2 + 1) / 

1)) + (^«) 



with = for i < i. 
It follows from (10) that 



For i — z ^ 0, this solution can be approximated as 

z 

^t^^) « (t) - 1, (12) 



as shown in Appendix B. 



2.2 Results of simulations 



We have conducted extensive simulations to evaluate the results presented 
in Section 2.1 on the in-degree distributions when the digraph evolves in a 
nearly acyclic fashion. The results we present in this section are averages of 
the quantities of interest over a large number of repetitions. 

Figures 1 and 2 refer to the case in which edges are deployed uniformly, that 
is, the Poisson-distributed number of edges that outgo from each node as it 
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Fig. 1. Average in-degree distribution in the nearly acyclic case for i = 10^ and 
z = 1,4,8, 12 when edges are deployed uniformly (10^ simulation runs). Solid plots 
give the analytic prediction of (5). 



enters the digraph is directed uniformly towards the nodes already present. 
Figure 1 contains plots of Pt{i, k) from both the simulations and the analytic 
prediction of (5) for different values of z and t, always with i — 10^. As the 
plots indicate, agreement is very good throughout. 

The plots of Figure 2 all refer to Pt{k), for which our analytic prediction is the 
one in (7). We show data for t — 10^ and several values of z. Notice that once 
again wc obtain good agreement between simulation and analytic prediction, 
but the two start to separate as k increases. The reason for this is clear: as k 
increases, the peak of the integrand of (6), occurring at x = kz/{z + 1), shifts 
continually to the right, which renders the extension of the integral's upper 
limit from the z\nt of (6) to infinity ever less justifiable. 

Data for the case of preferential attachment, in which edges get directed to- 
wards nodes already in the digraph with probabilities that arc proportional 
to how many incoming edges those nodes already have, are given in Figures 3 
and 4. Figure 3 is given for t — 10^ and a set of different z values, and confirms 
our expectation that in this case Pt{k) is a power law. 

For i = 10'^ and a few different values of z, Figure 4 contrasts simulation data 
for z^{i) with the analytic prediction of (12). Evidently, the two agree very 
well. 
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Fig. 2. Average in-degree distribution in the nearly acyclic case for t = 10^ when 
edges are deployed uniformly (10^ simulation runs). Solid plots give the analytic 
prediction of (7). 
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Fig. 3. Average in-degree distribution in the nearly acyclic case for t = 10^ when 
edges are deployed proportionally to the nodes' in-degrees (10^ simulation runs). 



3 Evolution allowing for cycles 

3.1 Analytic results 

Our second evolution scenario allows for the appearance of directed cycles 
other than self-loops by incorporating a probability, denoted by Pr, to control 

the replacement of certain edges by certain others. At each time step t > 0, 
one of the following two actions is selected to take place, recalling that the 
initial number of nodes is no = and that, as before, nodes are numbered 
consecutively from 1 as they enter the digraph: 
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= z = 4 




t 



Fig. 4. Average in-degree in the nearly acyclic case for i = 10^ and z = 1,4, 8, 12 
when edges are deployed proportionally to the nodes' in-degrees (10^ simulation 
runs). Solid plots give the analytic prediction of (12). 



• With probability pr, the number of nodes in the digraph remains the same, 
thus rit = rit-i- In addition, a node i for which df_i{i) > is selected and 
one of its incoming edges, chosen randomly with a uniform distribution, is 
replaced by a new edge directed away from i. Letting n^Li be the number of 
nodes that have nonzero in-dcgrccs at time t — 1, the probability that node 
i is selected is l/nf_i, provided its in-degree is nonzero. 

• With probability 1 — Pr, a new node, say i, is added to the digraph, so 
rit — rit-i + 1. Then a new edge is deployed from i to each of the digraph's 
rit nodes with probability mm{z/nt, 1}. 

Readily, for Pr = this evolution scenario is the one we considered first in 
Section 2. For p^ > 0, what it does is to induce the appearance of directed 
cycles, and consequently of nontrivial strong components. 

Under these new rules, rit is no longer deterministically equal to t, nor is node 
i necessarily the node added to the digraph at time t = i. Instead, rif is now 
a random variable with expected value (1 —Pr)t, so the expected time step at 
which node i is added to the digraph is t = — Pr), even though this can 
happen as early as time t — i. Once node i enters the digraph and receives its 
Poisson-distributed number of outgoing edges, the process of edge replacement 
acts to alter this number of edges randomly, so we still expect the out-degree 
(i) to be Poisson-distributed with mean (i) = z. 

Analyzing the distribution of the in-degree df{i), though, is significantly more 
complicated and is achieved by setting up a system of finite-difference equa- 
tions to describe its behavior. As in Section 2, Pt(i, k) continues to denote the 
distribution of in-degrees for node i at time t with t > i and /c > 0. In what 
follows, we assume that node i is already present in the digraph at time t — 1. 
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We do this for the sake of simphcity, and indicate how to compensate for it 
when we return to this issue later in this section. 



We start with the k > case. With probabihty 1 — Pr, Pt{i, k) is the sum of 
two probabihties, each corresponding to one of the following mutually exclusive 
events: (i) the node added to the digraph at time t connects out to node i; (ii) 
the node added to the digraph at time t does not connect out to node i. With 
probability Pr, three other mutually exclusive events must be considered: (iii) 
node i is the node selected to have one of its incoming edges replaced; (iv) 
node i is not the node selected to have one of its incoming edges replaced, nor 
is the replacing edge incoming to i; (v) node i is not the node selected to have 
one of its incoming edges replaced, but the replacing edge is incoming to i. 

Let us assume that, when node i is the node selected to have one of its incoming 
edges replaced, the probability that the randomly chosen replacing edge forms 
a self-loop, or connects out to a node towards which an edge from i already 
exists, is negligible. This is certain to hold as the digraph acquires more nodes, 
and allows us to conclude the following. In case of event (i) or (v), the in-degree 
of node i is increased by 1 from time t — 1 to time t, while in cases (ii) and (iv) 
it remains the same and in case (iii) it is decreased by 1. For k > 0, Pt{i,k) is 
then such that 



where we have used the fact that rit = Ut-i in the case of edge replacement at 
time t. For ntn^_i S> 1, and approximating rit by its expected value (1 — Pr)t, 
(13) becomes 



Pt{i, k)^{l- Pr) [mm{z/nt, l}Pt^i{i, k - 1) 

+ {l-mm{z/nt,l})Pt-i{i,k)] 




(13) 



Pt{i,k)-Pt-i{i,k) 




{l-Pr)Ht + 



Pt-l{i,k+l) 



Pr 



+ 



{1-Pr)t 



Pr 



] 




with /It — min {z/{l — Pr)t, 1}. 
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For the case of A; = 0, the only meaningful possibihties are events (ii)-(iv), 
hence 



n. 



t-i/ 

-((i-f>'+(r3^)^'.-.ft<')- (15) 



Note in addition that, for t > 0, it is also possible to approximate nf by its 
expected value, that is, 

nt ^nt-fl PtU, 0)^(1- Pr)t - E PtU, 0), (16) 
j=i j=i 



where the summations give the expected number of nodes that have zero in- 
degree at time t. The system of finite-difference equations given for /c > 
by (14)-(16) can then be solved numerically for alH > 1 and all i < t. The 
greatest possible in-degree for fixed i and t is t — i + 1, corresponding to the 
case in which node i is added to the digraph at time step i and acquires a 
new incoming edge at every subsequent step, so /c G {0, . . . ,t — i + 1}. All 
boundary conditions can be set to zero, and the value of Hq is immaterial, so 
long as it is set to some nonzero constant. 

In order to compensate for the fact that the finite-difference equations being 
solved for time t and node i are based on the assumption that node i is in the 
digraph at time t — 1, and also to ensure that a nontrivial solution is obtained, 
we may heuristically correct the equations for all t > 1, all i < t, and every 
appropriate k, by adding adequate probability terms to reflect the appearance 
of edges. If iitii) is the probability that node i is added to the digraph at time 
t, then we add 5ifc7rt(i) min{2;/i, 1} to (14), where 6ik is the Kronecker delta 
function for k = 1, and T^tii) (1 — mm{z/i, 1}) to (15). 

As for 7rt(i), it is given by 



and thus admits an approximation based on the Poisson distribution, that is, 
7rt{{) ^ (l-p.)P((l-p.)(t-l),^-l), (18) 

where we recall that P{x, k) denotes the Poisson distribution with mean x (cf. 
[1])- 
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Fig. 5. Average out-degree distribution under directed-edge replacements for 
i = 10^, z = 1,4,8, 12, and Pr = 0.25,0.75 (10^ simulation runs). Solid plots show 
the Poisson distribution with mean given by the corresponding value of z. 



3. 2 Results of simulations 

We start with a presentation of simulation data intended to be confronted 
with our predictions in Section 3.1 regarding the distributions of in- and out- 
degrees when, during its evolution, the digraph is allowed to have directed 
cycles other than self-loops by the action of edge replacements. These data 
are presented in Figures 5 and 6. 

Figure 5 contains plots for pr = 0.25 and pr = 0.75, along with a few values of 
z and t. Each plot shows simulation data on the distribution of the out-degree 
df{i) of node i — 10^. In the figure, we use Qt{i,k) to denote the proba- 
bility that dt{i) — k, and also show the corresponding Poisson distribution, 
P{z,k). As we mentioned earlier, our expectation is for Poisson-distributed 
out-degrees, because the number of outgoing edges a node is given upon en- 
tering the digraph is thus distributed, and from there onward all that the 
edge-replacement mechanism does is to randomly alter the initial deployment 
of edges. We see in Figure 5 that this is indeed the case to a very good degree 
of agreement between Qt{i, k) and the Poisson distribution with mean z. 

Data for the in-degree distribution Pt{i, k) are shown in Figure 6 for i = 10^ 
and two values of {Pr — 0.25 andp^ — 0.75). All plots are shown as functions 
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Fig. 6. Average in-dcgree distribution under directed-edge replacements for i = 10^, 
z = 1,4,8,12, and p,. = 0.25,0.75 (10^ simulation runs). Panels placed side by 
side share the same value of z. Solid plots give the solution of the system of fi- 
nite-difference equations in (14)-(16). 

of t for a few different values of z and k. For each combination of z and k 
values, two plots are given, one for simulation data and another to depict the 
solution of the system of finite-difference equations in (14)-(16) that for A; > 
describes the behavior of Pt{i, k) for fixed Pr and z. It is clear from the figure 
that agreement is very good between the two plots in all cases, with the very 
few exceptions of some of the higher values of k, but such discrepancies are 
the result of insufficient statistics for those particular values of k, despite the 
10^ simulation runs. Agreement is so good, in fact, that solving the system of 
equations becomes largely preferable to simulating the system a sufficiently 
large number of times, since the latter is unavoidably slower by several orders 
of magnitude. 

A by-product of our simulations has been the empirical finding that the ex- 
pected vahie of n^" is [z/{z + 1)](1 —Pr)t, where we recall that (1 — Pr)t is the 
expected value of n^. This fact is very useful, since it decouples the calculation 
of n^" from that of the in-degree distributions of all nodes (cf. (16)). Conse- 
quently, it becomes possible to solve the equations in (14) and (15) for only a 
selected set of nodes of interest. 

We now present a series of four figures. Figures 7-10, where we show data 
related to the behavior of the digraph's strong components, which under the 
policy of edge replacement exist in nontrivial form. All the plots given in 
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Fig. 7. Average number of strong components (NSC), of cycle components (NCC), 
and of knots (NK) as a function of z for pr = 0.25,0.75 (10^ simulation runs). 



these figures are shown as functions of z for p,. = 0.25, Pr = 0.75, and a few 
relevant t values. The data presented henceforth are then to be contrasted 
with the expectation we have from the static case studied in [1], since the 
dynamic case is expected to become qualitatively equivalent to the static case 
for large enough t, so long as > (the larger the value of pr, the sooner the 
equivalence can be observed as t grows). The intuition behind this statement 
is that, when edges are continually replaced, the passage of enough time steps 
is expected to add such a degree of randomization to the positioning of the 
edges that it becomes probabilistically indistinguishable from any positioning 
that could result if edges were deployed totally randomly to begin with, as in 
the static case. 

Figure 7 is devoted to showing how the number of strong components be- 
haved during the simulations we conducted, and also in particular the number 
of cycle components and of knots. In the static case, the number of strong 
components tends to one quickly right past z = 1 as z is increased and the 
giant strong component appears. In the dynamic case we expect a similar ef- 
fect as both z and t are increased, that is, we expect the number of strong 
components to be dramatically reduced. This is what the figure demonstrates, 
particularly for the higher value of Pr- 

Regarding the number of cycle components, in the static case it peaks at 
z — 1 and decreases rapidly to either side, eventually approaching zero. The 
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data in Figure 7 tend to support the expected peak at z = 1, but even for the 
higher vahic of Pr a nonzero number of cycle components, albeit small, seems 
to be sustained as both z and t grow. We will see shortly that this is due to 
the presence of cycle components that are self-loops, and presume that it can 
be explained by demonstrating that such components have a non-negligible 
probability of appearing for small values of t, being on the other hand less 
likely to be picked for replacement as t increases. 

In the static case, knots are expected to comprise one single node and to occur 
sparsely for very small values of z, then to be altogether absent as z is increased 
but still kept below ln(2n), and then to occur as a single, all-encompassing knot 
for z > ln(2n). In the dynamic case, wc expect a similar behavior, including the 
effect around the threshold z = ln(2nt) ~ In (2(1 —pr)t). Taking pr = 0.75 
and t — 8000, for example, yields z — 8.3, which agrees well with what is 
shown in Figure 7. But the absence of knots for the intermediate values of z 
seems too brief when contrasted with the expectation created by the static 
case. Returning to the corresponding data in [1], specifically Figure 6, we see 
that our expectation was built on one single simulation run, which may have 
been the cause for the seemingly larger intervals of z inside which knots were 
totally absent. 

These observations are complemented by the data shown in Figures 8-10, 
where minimum and maximum sizes are shown for the strong components in 
general, the cycle components, and the knots, respectively, observed during the 
simulations. As shown in Figure 8, the minimum size of a strong component is 
one, particularly as t grows, while the maximum approaches Ut ~ (1 — Pr)t as 
z grows away from one, reflecting the expected appearance of the giant strong 
component comprising all rit nodes. 

Minimum and maximum sizes observed for the cycle components are as shown 
in Figure 9. Aside for a small vicinity near z = 1, where sizes larger than one 
occur, cycle components have size one for all values of t. This lends support 
to our earlier suspicion that the cycle components still persisting as both z 
and t grow are in fact size-one components. 

The case of knots is illustrated by the data in Figure 10. It shows that the 
single knot that appears suddenly roughly around z — hi (2(1 — Pr)t) is indeed 
all-encompassing, as its size tends to rit ~ (1 —pr)t as z grows, thus confirming 
our expectations from the static case. 
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Fig. 8. Average minimum and maximum sizes of strong components (SC) as a 
function of z for pr = 0.25,0.75 (10^ simulation runs). 
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Fig. 9. Average minimum and maximum sizes of cycle components (CC) as a func- 
tion of z for Pr = 0.25, 0.75 (10^ simulation runs). 
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Fig. 10. Average minimum and maximum sizes of knots as a function of z for 
Pr = 0.25,0.75 (10^ simulation runs). 



4 Conclusions 



We have in this paper considered random digraphs that grow in discrete time 
by the continual addition of new nodes and edges. Within this context, we 
started with a study of digraphs that remain nearly acychc at all times as 
nodes are given a Poisson-distributed number of outgoing edges upon entering 
the digraph. For two modes of attachment to the nodes already in the digraph, 
uniform and preferential, we contributed analytic and simulation results that 
describe the digraph's degree distributions. 

In the second part of the paper, we allowed edges to be randomly replaced 
during the evolution of the digraph, aiming at allowing nontrivial strong com- 
ponents to appear. In this case, too, we contributed analytic and simulation 
results, with emphasis on the introduction of a system of finite-difference equa- 
tions for the computation of a node's in-degree distribution at all times. We 
finalized by returning to the strong components that were our main subject 
in [1] and investigated their appearance and behavior along the evolution of 
the digraph. 

We mention that it is possible to extend our study of the evolution of a digraph 
and its strong components to, in principle, all kinds of connection rules and 
evolution strategies, particularly those that have special significance within a 
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certain application area. The ones we studied were ultimately targeted at the 
study of strong components, and for this reason the absence or presence of 
directed cycles, as well as the conditions for their appearance during the evo- 
lution, have been crucially important. Other guiding principles will certainly 
exist as the focus is shifted by some other motivation. 



A Derivation of (3) 

Note first that, for each subset of {i, . . . ,t} comprising k nodes, Pt{i, k) ex- 
presses the probabihty that an edge exists directed from each of those k nodes 
to node i, but not from any of the remaining t — i + 1 — k nodes of {i, . . . , t}. 
Let C denote the set of all partitions (C, C) of {i, . . . , t} such that C contains 
exactly k nodes. Then we have 




ic,c)ec xec 





(A.l) 



Since Pt{i, 0) 



n^=i(l — z/x), we can write (A.l) as 




(A.2) 



For -C t — i, in (A.2) we can use the approximations 




(A.3) 
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and 



E n 1 

BCCUC X^B ^ 
\B\=j 



X- ZJ J 



xecuc ^ 



X — z 



t-i + l + ln(- 



t - z 



\i — z — 1 



(A.4) 



which yield 



k\ 
k\ 

and, for z <^ i,t, 



j=0 



In 



t-i + l + ln( 



t - z 
\i — z — 1 



t- z 
i — z — 1 



ln( T 



z-i k 



(A.5) 



(A.6) 



by (2). 



B Derivation of (12) 



From (11), we have 



x=i \ 



x{z + 1) 



- 1 



t-i+1 



E 



1 



j=o + i<xi<-<xj<t ^1 ■ ■ ■ 



1. 



(B.l) 



For i — i ^ 0, in (B.l) the innermost summation can be approximated by 

t 



In 



i-1 



(B.2) 



thus yielding, for i ^ 0, 



~ E 7T 



t-i+1 



ln(f^ 



(B.3) 
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